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Dark matter annihilation and other processes may precipitate a flux of diffuse ultra-high energy 
7-rays. These 7-rays may be observable in present day experiments which observe diffuse fluxes 
at the GeV scale. Yet the universe is presently opaque to 7-rays above 10 TeV. It is generally 
assumed that cascade radiation is observable at all high energies, however the disparity in energy 
from production to observation has important consequences for theoretical flux limits. We detail 
the physics of cascade radiation development and consider the influence of energy and redshift scale 
on arbitrary flux limits that result from electromagnetic cascade. 



I. INTRODUCTION 
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FIG. 1: The pair-photon cascade cycles 7-rays through a series of steps that cool these 7-rays while conserving energy overall. 
The total injected power is observable in output spectra. First, diffuse isotropic sources inject 7-rays above the threshold energy 
(1). Next these 7-rays cool through pair production and inverse Compton scattering (2). Finally, the cooled 7-ray spectrum is 
observable as a spectrum of inverse Compton photons (3). 



The GZK process produces a flux of ultra high-energy 7-rays. These gamma rays are scattered by the CMB. 
Furthermore, recent work concerning dark matter annihilation into 7-rays suggests this is an important pathway for 
constraining and revealing the nature of dark matter particles [l| . Evidently the universe is replete with non standard 
model particles. These particles must account for roughly 90% of the matter budget for the universe as a whole. 
While the galaxy scale distributions of these particles remain in contention, it is known that on large scales dark 
matter particles are distributed with relative uniformity. Any uniform distribution of particles that annihilate to 
7-rays are subject to substantial observational constraints in both existing and planned experiments. A caveat occurs 
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when these particles are produced above the threshold energy for pair production at cosmological scales. The flux of 
particles at Earth may be mitigated by cosmogenic propagation. Therefore, before one may constrain particle fluxes 
based on electromagnetic cascades, one must have a detailed understanding of cascade development, energy scale, 
production time scale and detector sensitivity. 

In addition, the origin of the diffuse extragalactic component of EGRET observations is a deeply held mystery in 
7-ray astrophysics @, y|, E| • In the galaxy it is predicted that inverse Compton scattering of cosmic-ray electrons on 
interstellar photons may be the source of diffuse 7-ray emission || . Alternatively, both cosmological and exotic 
ultra-high energy processes may contribute a portion of the extragalactic 7-ray spectrum. Fermi shock acceleration 
may inject ultra-high energy cosmic-ray electrons beyond the TeV scale All of these extragalactic high-energy 
processes are rare and diffuse. If these processes contribute 7-rays above the pair production threshold, then the 
determination of the observable spectra resulting from diffuse isotropic cosmological injection is an essential precursor 
to considering astrophysical models. 

II. THE CASCADE 

An isotropic injection of ultra-high energy 7-rays is subject to the pair-photon cascade. The cascade was suggested 
as early as 1948 [9j, however it was described in detail by Bonometto and Rees in 1971 [10]. This process results when 
the coupled chain reaction of pair production (JTJ) and inverse Compton scattering @ is possible. 

77 — * e + e~ (1) 



e7 — > ej (2) 

Berezinsky illustrates the basic argument in several works [ill . fl2L Il3j |. First, an incoming high energy 7-ray 
encounters a thermal background photon and forms an e + e~ pair. Pair production is the subject of a useful review 
by Motz, Olsen, and Koch [1J|. At low energies, with s-wave scattering, the outgoing electron and positron share the 
incoming total energy: 

E e = E 7 /2 

At high energies, a leading particle carries away most of the incoming energy. This is the forward scattering limit. 
The threshold for pair production presents an absolute lower bound on the energy where the cascade may occur. This 
threshold depends on the average energy of isotropic target photons, eo, and the mass of the electron. At present this 
value is: 

E th =ml/e Q = 10 (^) TeV 

At this energy the target distribution is given by the infrared background, while at PeV leading energies the cosmic 
microwave background may participate. For linear cascades, which develop in low density regions of the universe, the 
products of cascade development do not participate directly, therefore it is safe to assume that the target distribution 
is mono-energetic, as in the following brief discussion. Later the thermal photon spectra is presented along with the 
kinetic solution. Stawarz and Kirk discuss non- linear cascades in a recent work [151 ]. 

In all of the past discussions of cascade development, one crucial factor is typically under represented. The back- 
ground distributions of target photons typically evolve rapidly with redshift, e = eo (1 + z) [16]. We will see that 
redshift plays an important role in observability of cascade fluxes. 

After pair-production occurs on target photons, the resulting stream of high energy electrons and positrons are 
susceptible to inverse Compton scattering. Through inverse Compton scattering the energy of these electrons or 
positrons is transferred into a background photon which results in a new high energy outgoing 7-ray. The electrons 
and positrons are left behind and do not participate further. The vast majority of the incoming energy is carried by 
the outgoing 7-ray. In the low energy limit, the average energy loss fraction is [i~7l ]. 
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An outgoing electron or positron resulting from 7-ray injection at threshold will have energy: 

Ee = ^m 2 Je 

Thus, the Cascade will exhibit a transition energy when outgoing 7-rays no longer have sufficient energy to initiate 
pair production. This transition energy, called the critical energy, E c , is the energy of an outgoing inverse Compton 
scattered 7-ray resulting from an electron or positron produced with energy Efh/2. 

or 



E c — E ltOUt — -E 2 e /ml — -E 7 j h 



At present, the cascade begins to transition from recycling to emission above 3 TeV. This process is depicted in 
Fig. [TJ The net effect of the cascade is the reprocessing of one incoming high energy 7-ray into pairs of outgoing 
7-rays, each with roughly half of their parent's energy. These pairs form others until numerous final particles result. 
All of these final particles combined share the energy of the incoming 7-ray. 

There are two important energy scales for cascading particles. Below the energy of transition and above. Above E C1 
energy is conserved by the cascade process, therefore we may immediately write the spectrum, E 2 dN/dE — const. 
Below E c , the last generation of 7-rays are unable to pair produce, therefore they must escape after inverse Compton 
scattering which has no low energy threshold. These 7-rays lose energy and conserve number. A deduction of this 
spectrum is given in App. fA] 



(4) 

This conventional approximation to the cascade is prevalent in many studies of 7-ray fluxes. The analysis presented 
considers the effect of one step of the cascade and allows for reasonable approximations to be made for given models. 
It is also very simple to consider the effects of two cascade steps. This derivation is given in App. [BJ Finally, for a 
multi-stage cascade the following approximation is appropriate: 

,,. \AE- 2 E>E th 

^1 = {A'E-^ E c < E < E t h (5) 
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The cascade conserves total energy as it cools the incoming 7-rays and one exploits this to fix the relation between 
injection spectra and emission spectra and to determine the normalization constants. 

This discussion is essentially complete less one crucial detail. While present observations exist at GeV energies, 
the present value of E c is about 10 TeV. In order to develop limits on particle fluxes we must consider this issue. 
The opacity of the universe evolves in redshift. Present experiments will have sensitivity to electromagnetic cascades 
which develop at z 1 and beyond. If cascade radiation is produced only at present, it may only connect with present 
diffuse observations through extreme downscattering. 

All cascade flux limits are deduced from the assumption that the energy density of cascade radiation does not 
exceed experimental sensitivity as in Fig. [2j For an energy density p [GeV/cm 3 ], the following relation is true today: 

Peas < Pobs (6) 

The flux limit on particles injected above E c directly follows from this relation: 
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FIG. 2: EGRET is sensitive to extragalactic 7-ray backgrounds below 30 GeV (left) 4]. When compared with fluxes of 
cascading particles it is evident that experimental limits made at present apply only above 2fsl (right). 

The most important relevant observation is from EGRET, see Fig.O The EGRET experiment on-board the Compton 
7-ray observatory has set an upper bound at 30 GeV of roughly 10~ 6 GeV cm -2 s _1 sr _1 Q. 



Therefore there are two ways to develop cascade limits based on observations at present. Either cascade limits 
must be deduced from a prior epoch when the universe was opaque to 7-rays at 30 GeV, i.e. z w 1, or limits must 
incorporate the emission spectrum of 7-rays below E c to decide what part of the present era flux contributes to a 
particular experiment. 

First, it is straight-forward to predict the redshift effects on cascade radiation. Consider, the cascade radiation in a 
comoving volume V defined for an arbitrary but fixed number of cascade particles. Then by ([5]), p cas V < PobsV is also 
true. However, V may be parameterized by a scale factor a(t), and coordinate representation F(r), V = A/ina 3 F 3 . 
Since a parameterizes all time dependence and similarly F all spatial dependence one may easily conclude that, 
Pcasa 3 < PobsO,Q is always true. Here ao is used to reference the present day value of the scale factor. This is more 
customarily written in terms of redshift: 



^ I a \3 Peas 
Pobs > Pcas{ ) ' 
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These redshift effects are generally canceled by source evolution with density, i.e. n = no(l + z) 3 . 

Alternatively, cascade radiation may be contributed by downscattering of final stage cascade electrons and positrons. 
We will treat this spectral development in detail below, yet one can make an initial estimate by taking inverse Compton 
scattering to have an E~ 3 / 2 spectrum. By (|6]), for emission below E c and observation at E {, s one has p b s N 7 > p cas N 7 



t-,-3/2 
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> PcasE, 
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As in the case of EGRET observations, if E c /E Q t, s w 100, then the correction factor will 



be 1/1000. Even GLAST will only reduce this to 1/100. Unfortunately, present day cascade radiation will almost 
certainly not contribute to present day experimental observation through this mechanism. Eq. ([6]) must be extended 
with the following relation: 



PcasiE) = 



dE' I dz 
'0 



Pcas(E c ( Z ) = E') (E, 

(1 + z) 3 \E< 



-3/2 



(7) 



From this one may directly consider a few simple cases: 



1. Present era high energy injection and observation at E c : In this case p' cas = p ca s(E c ). 
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2. Present day high energy injection and observation somewhat below E c . In this case one has a (E c /E b s ) 3 / 2 
penalty: 




e c \- 3 / 2 



For the currently practical case of EGRET, contributions to cascade limits by these fluxes are ruled out by a 
factor of 1/1000. This could only result in very weak limits on particle fluxes. 

3. Observation of a 30 GeV process with significant contributions at z = 1. This might be considered the 'standard' 
or conventional case. Here the best case contribution is: 



Consideration of other processes requires a more detailed discussion, however it is already evident that 7-ray cooling 
has important consequences for cosmological production processes. Any 7-rays injected above E t h may be observable 
below this energy. The source of the diffuse flux of 7-rays below 30 GeV is an open question for astrophysics. However, 
The energy resulting from cascade 7-rays is narrowly distributed around E c with little or no contribution as energies 
approach 30 GeV. These fluxes are suppressed by three orders of magnitude. Important cosmological processes would 
have to inject enormous amounts of energy above E t h to even partially contribute to EGRET diffuse observations 
today. Therefore it is unlikely that the source of the EGRET diffuse observation is local 7-ray injection above E t h- 



Let us now turn to the saturated pair-photon cascade problem to develop a more detailed spectral understanding of 
cascading particles. For cosmogenic processes one is only concerned with saturated propagation. High energy injection 
processes are defined to be saturated if through extremes of cosmological scale or density all injected 7-rays eventually 
must pair produce 20]. In short, the universe is completely opaque to these particles. Obviously, this calculation is 
not sensitive to conditions in which the universe is not completely opaque. The cascade does not conserve particle 
number so repeated solution of scattering relations is required to determine a final emission spectrum. 

To solve saturated propagation it is customary to employ the method proposed by Guilbert [lq . An integral kinetic 
equation gives the spectra resulting from sources. The cross sections of 7-ray scattering are well known. It is possible to 
directly calculate the probability of a 7-ray scattering from a given energy to any other energy. Repeated convolution 
of the scattering probability with the incident number distribution finally determines the resulting spectrum. 

Svensson and Zdziarski [19|, [2(| HH HH, [23|, [24], HH, [26[ considered the cascade problem extensively using a similar 
analytical approach. The steady-state solution to a kinetic equation of propagation allows an analytical determination 
of final spectra based solely on rate and injection models. If electron escape is neglected, then the steady state electron 
distribution is completely described by electron production and energy loss. 

This solution applicable to diffuse cosmological processes makes the following assumptions. 

1. A narrow isotropic distribution of 7-rays is injected above E t h- This may similarly be solved for power laws and 
other distributions. 

2. The interaction scale for photons and electrons is short (10 Mpc) in comparison with propagation scales (Gpc). 
Therefore, cosmogenic injection processes are saturated above the pair-production threshold, E t h- 

3. The universe is pervaded by an isotropic distribution of soft background photons at a recent epoch, z < 0.03. 
This solution may be extended to higher redshifts without difficulty. 

4. The energy loss time-scale for electrons is negligible in comparison with any possible escape time-scale. 

5. Any homogeneous magnetic fields present on cosmological scales have a negligible effect on energy loss. 

Based on these assumptions, a self-consistent solution to the saturated pair cascade problem in the ultra-high energy 
low redshift regime follows. The observable 7-ray spectrum below E c is then deduced. It turns out that the isotropic 
cosmological background acts as a photon calorimeter, while the conformation of any injected spectra are lost, total 
energy of injection is preserved and observable in experiments with sensitivity near E c . 



Peas 




III. SATURATED PAIR CASCADES AND DIFFERENTIAL FLUXES 
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IV. INTERACTIONS OF COSMOLOGICAL 7-RAYS AND ELECTRONS 

Cosmological 7-rays and electrons may be susceptible to energy loss through inverse Compton scattering, pair 
production, photon-photon scattering, Compton scattering, synchrotron radiation and redshift. One may briefly 
summarize the effects of these processes to deduce the important propagation processes. In the conventional notation 
for dimensionless energy, one refers to photon energy with w — E 7 /m e and electron Lorentz factor 7 = E e /m e . 



Incident Particle 




Process 


Target Density [cm ' 


'] Cross Section [<jt] 


A [Mpc] 


Gamma Rays 
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v e- on CMB (PeV) 
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3/16 
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Electrons and Positrons 




e7 — ► 


ey on CIB (TeV) 


0.5 


1 


1 




e7 -» e7 on CMB (PeV) 


410.5 


1 


10~ 2 



TABLE I: Inverse Compton scattering and pair production are important cosmological loss processes for electrons and 7-rays. 



1. Pair Production from y-Rays 

To estimate the scattering length for 7-rays I assume cosmological 7-rays will pair-produce on an isotropic infrared 
background if e 7 w 10 TeV. The center of mass energy squared for TeV gamma-rays against infrared background 
photons is s w 4e 7l e 72 ~ 10~ 13 TeV 2 , the dimensionless center of mass energy is ^/UJilo^ ~ 1. The cross section is 
found at the peak of the pair production rate a w jq&t- The Thomson cross section, <tt, is related to the classical 
electron radius by 87rr 2 /3. The mean interaction length for cosmogenic pair production is about 10 Mpc. This mean 
free path is minimized for injected gamma-rays at PeV energies. In this case, the interaction length drops to about 
10 kpc because of the increased density of microwave background targets. 

A~~_ > .+.- = — w 10 Mpc 
77 no 

In comparison with the Hubble scale, pair production is a primary energy loss mechanism for cosmological 7-rays 
with energy above E t h- 

2. Compton Scattering of Gamma-Rays 

Compton scattering of a cosmological 7-ray on a primordial electron is primarily important at epochs where primor- 
dial electron density is significant. A TeV 7-ray and background electron system has center of mass energy squared, 
s w 4e 7 m e w 10~ 6 TeV 2 . In dimensionless form, y/uJy rj 3000, the cross section for TeV 7-rays is a ~ 10~ 3 ar~ 
10~ 5 <tt- The density of primordial electrons is roughly equivalent to the density of baryons, n e w rib ~ 10 -7 cm -3 . 
The mean interaction length of Compton scattering is only significant at high redshifts. 



A 7e ^ 7e = — rj 10 10 Mpc 
no 

The mean free path for Compton scattering is so large that I take these losses to be negligible. 
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3. Photon-Photon Scattering of Gamma-Rays 



Photon-photon scattering, 77 — ► 77, is an important cosmological consideration at high redshift and low energy 
scales when the universe was radiation dominated. Svensson and Zdziarski treat this process in detail in [25] . I neglect 
this scattering process here under the assumption of ultra-high energy propagation in a matter dominated epoch. 



4- Inverse Compton Scattering of Electrons and Positrons 

Cosmological electrons of TeV energies may scatter on isotropic primordial photons. The center of mass energy 
squared is s ~ 4£' e e 7 . The electron energy, E e , is about 1 TeV. The energy of infrared photons averages about 
10 -2 eV, so the dimensionless center of mass energy is « 0.5. This is well within the Thomson regime so the 
cross section is the Thomson cross section. The density of soft infrared target photons is 0.5 cm -3 . The interaction 
length of TeV electrons is about one Mpc. This is a rough estimate made for illustrative purposes, the integration 
below uses exact distributions. This interaction length is minimized when the CMB can participate. At PeV energies 
the interaction length of an electron drops to about 1 kpc. 

A e7 ^e 7 = ~ 1 MpC 

tiut 

As a result inverse Compton scattering is a primary energy loss mechanism for cosmological electrons above one TeV. 



5. Synchrotron Radiation from Electrons and Positrons 



Cosmological electrons in a magnetic field, B, emit 7-rays as they accelerate through a helical trajectory. The 
energy loss through this mechanism is proportional to magnetic energy density or B 2 . 



dE 4 2 2 B 2 

The cosmological magnetic field density has been constrained at less than 10~ 12 Gauss. Since synchrotron loss depends 
on B 2 , 1 assume that cosmological magnetic fields may be neglected in this calculation. See [27J for a detailed treatment 
of this loss process. 



6. Redshift 



The energy loss scale due to redshift is of roughly the same order as the Hubble length scale. 



A, « 4" ~4G P c 
no 

The extension of this solution to include detailed accounting of redshift is straightforward but not required since the 
principal loss processes are far faster than this mechanism. 



V. PRODUCTION IN ASTROPHYSICAL PROCESSES 



For particles susceptible to these interactions a particular isotropic injection has flux defined by 



$ = d E— (8) 

Air J Emin dE 1 ' 

The differential number density is taken as dN/dE [GeV -1 cm -3 ]. Where dN/dE dE represents the number of 
particles in an interval [E, E + dE) per volume. If N represents the density today, the density evolves with redshift, 
z. 
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N = N (l + zf 

The flux of particles resulting from a particular astrophysical process is determined by the production rate. Rate 
is given by ((§]), expressed in terms of the incident particle distribution n, and reaction cross section, a. 

T(E',E)=c den(e) dE'—(E',E) (9) 

J (min JE min dE 

The term ^(E',E) represents the differential cross section for scattering from a particle of energy E' to energy 
E. Also this integral can be normalized over target particle distributions and arbitrary directions. This paper adopts 
the notation \x = cos for the scattering angle and uses j3 to refer to velocity. The following equation represents the 
production from scattering of an isotropic incident distribution on a target distribution across a spectrum of energies. 



•V = r I ,h n ri<) i <!fi I dE! nj(Bj) f """ dE' ^(E',E) (10) 



Ei 



In this equation ut represents the distribution of targets and nj represents the distribution of incident particles. I 
use N to represent the first derivative in time, dN/dt, hereafter. 



VI. THE KINETIC EQUATION 

The kinetic equation gives the observable spectrum, N(E), of particles after scattering for a given injection spectrum 
Nq(E). The transport of high energy 7-rays through the cosmological medium involves catastrophic energy loss 
through pair production. The resulting stream of electrons boosts photons from the primordial background through 
inverse Compton scattering. 

The contributions to 7-ray energy loss are either continuous or discrete. The continuous radiative transfer of a 
propagating particle is described as a differential equation over spatial propagation (34| . 



(E) = -aN(E)+j (11) 

dx 

Where a is the coefficient of assumed continuous energy loss and j is the particle injection term. Solutions to (fTTj) 
reveal the effects of continuous energy loss on an initial spectrum, but do not describe catastrophic loss processes. To 
consider both types of loss processes one employs a steady-state differential equation. 

The observable spectrum of 7-rays can be deduced by considering the repeated effects of energy loss in multiple 
scatterings. Since escape is neglected, differential changes in particle flux must be stable. 

The electron steady-state is described by loss due to inverse Compton scattering, production of e + e~ pairs, and 
injection (20j . 



Ne = ~N e ,c{E) + N e ,p{E) + N e ,in{ E ) (12) 

Likewise 7-rays appear after inverse Compton scattering or injection, and disappear in pair production. I will use 
n (lower case) to label a spectrum of photons and N to label a spectrum of electrons. 



n 7 = h y C (E) - h lt p(E) - C ~^( E ) + h ltin (E) (13) 

In (TT2")) and (fT3")) the subscripts "C" and "P" are used to indicate the time derivative due to inverse Compton 
scattering ^ and pair production ([T]) respectively. The radiative transfer is given by dn/dx and (fTTj) . The radiative 
transfer term may be used to model losses from 77 absorption in dense environments. The terms with "in" subscripts 
represent continuous isotropic particle injection. Both equations are independent of charge, therefore electron and 
positron losses are identical. In order to consider the effects of both electron and positron production one simply 
doubles the electron production rate. For the remainder of the article I will refer to electrons only, of course in reality 
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half of the denumerable electron population would be physical positrons. This is irrelevant since I only consider the 
resulting flux of 7-rays. 

For many physical models involving the injection of 7-rays there are no electron sources, in this numerical integration 
N e ,in — 0, this choice is arbitrary and made for numerical convenience, it would be equivalent to inject first generation 
electrons or positrons rather than 7-rays. 

Since (|12|) and (|13|) are coupled by pair production and inverse Compton scattering, one may adopt the point of 
view of either electrons or 7-rays, I choose electrons in consideration of previous work in this field. 

Now by considering the known cross sections for inverse Compton scattering and pair production one may write a 
differential equation describing electron production in terms of energy. Inverse Compton scattering on an isotropic 
distribution of background photons is given by (fT4")h 

N e ,c(E) = c /"""" de n T (e) [ dfi dE' N e (E')^-(E' ,E) (14) 



1 



If the background photons are thermal, the density is integrated from the Planck distribution. The present day 
temperature is To = Tqmb = 2.725 K and kb is the Boltzmann constant. In general, T is proportional to Tq(1 + z). 
The energy evolves with redshift according to the relation e = eo(l + z). 

n T (e) de = 1 -de 

ir 2 (hc) s exp(^) - 1 

The primarily relevant background for 7-rays at 10 TeV is the infrared background (CIB). I assume this background 
is adequately described by a power-law and a black body at Tcib = 2725K. These assumptions are consistent with 
recently published detailed models of photon backgrounds 1(| 28| . I arbitrarily normalize to achieve agreement with 



accepted observations, see citations above for detailed discussions of these backgrounds. 

The Compton energy loss rate of electrons is equal in magnitude to the energy production of Compton 7-rays. For 
a photon of energy ujra e c 2 , and an electron of energy r )m e c 2 , the exact angle averaged scattering rate for electron 
disappearance is given by the Klein-Nishina form [291 ]. 

iV e ,c(w,7) = c / den T (e) dn N e ( K )^^- (15) 



J2y(i—p)w dn 



-\ K ) 



dn v ' 32 7 2 /3u; 2 



4 8V, ,18 1 



1 2 Ml + «)+ o + -- 



k n 2 J y ' 2 K 2(1 + k) 2 



It is usual to employ a change of variables, k = 2j(l — /3/i)w, in this equation. In the low energy ( "Thomson" ) limit 
this rate approaches cax, the speed of light multiplied by the Thomson cross section. The Compton rate is depicted 
in Fig. H 



Equation (|15[) can be analytically integrated and this is a common step for several authors including Zdziarski, 
however the resulting rate is a complicated function involving dilogarithms. There may be no computational or 
intellectual benefit from performing this integration, I omit it. In either case one is required to numerically integrate 
over the rates to deduce a final spectrum. This approach is also suggested by Coppi and Blandford [29l ]. 

To clarify the meaning of (fTS")) I introduce a shorthand, C(E) = — dE'da/dE, for the portion of the integral 

over (TTJ} from E m i n up to the energy of consideration E. It is evident from Fig. [3] that the derivative of the cross 
section is negative. 

N e ,c(E) ^ c J^ ma * de nT(e) J dfi l^M Ne ( E ) 

x (c(E) - J^* dE' ^§{E',E)\ (16) 

In (jT5J) the first term represents electron energy loss due to boosted primordial photons, it is positive (loss) because 
these inverse Compton electrons are scattered to lower energies. The second term represents the appearance of inverse 
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FIG. 3: The scattering rate of 7-rays due to the inverse Compton process Electron energy is fixed at 20m e c 2 . This figure 
corresponds to Coppi and Blandford Fig. 1 [29|: given here for comparison. The Compton rate is unusual in that it decreases 
as energy increases. The maximum rate in the "Thomson" regime occurs at car- The rate is divided by car to reflect the 
absolute shape of the interaction in dimensionless form. 

Compton scattered electrons from higher energies. The sign is consistent with (1121) which defines this equation as a 
loss rate. 

The pair production or electron appearance rate is given by 



NeAE) 



den T (e)J dfi J dE y n 7 (Ey) -j-^(E 7 , E). 



(17) 



The term ^^-{E^^E) represents the differential cross section of a photon of energy E y to produce an electron of 
energy E. Eq. (flU)) reveals that there are two sources of 7-rays that may take part in pair production, either freshly 
injected 7-rays from isotropic sources or 7-rays boosted through inverse Compton scattering. 



n y (E) = n Jtin (E) + n y ^{E) 



These two relations are combined in (fT5)) . 



N e (E)= I den T (e) C dfi 



dE y ~T&( E i> E ) 
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(18) 



The exact photon-photon pair production rate is given by [2E 



1 — fEmax ^ a 
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FIG. 4: The scattering rate of a photon of energy io\m e <? against a photon of energy W2n.ec 2 to form electron-positron pairs 
|T|). This figure corresponds to Coppi and Blandford Fig. 5 [2y]: given here for comparison. The pair production rate peaks at 
about 3<tt/16 and then falls off as energies increase. The rate is divided by cctt to reflect the absolute shape of the interaction 
in dimensionless form. 

Following standard notation, dimensionless parameters for photon energy are u>i — e/m e c 2 and 0J2 = E 7 /m e c 2 , 
so [imax — max(— 1, 1 — 2/(jj\U2)- The electron velocity is {3' = [1 — 2/<JiW2(l — fi)] 1 ^ 2 . The pair production rate is 
depicted in Fig. |H 

For extragalactic cosmic media it is assumed that the inverse Compton scattering and pair production mechanisms 
entirely describe the cascade energy loss and no additional energy losses are present, dN/dx — 0. Combining (| 1 2[) . 
([in]) and l|18p reveals a steady-state integral-differential equation of electron transport. 



N e = -c 



de n T (e) J dp, } 



N e (E)(c(E)-J^dE' d -^(E',E) 



l-/x ' E * 



d< „,(,] j d„ j dE, ^§(E„E) 



X fl yi in(E) + c 



E C ,, 



Ec.r, 



dE' N e (E')^-(E',E' -Ey) 



(20) 



Eq. (|20[) is comparable to Zdziarski Equation 1 [21(, with solutions that show the time evolution of the electron 
spectrum versus energy. It is possible to eliminate time dependence by considering the continuous energy loss dEc / dt 



Here, Zdziarski gives the continuous energy loss in terms of a small parameter S, I take S to be 10 -4 [2l| : 



dE c 
dt 



(E) 



£7(1+5) 



dE' (E' - E)N e> 



c 
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This alteration gives an equation comparable to Zdziarski Eq. (2) [2lj : 

d ' dEc N e ]=-c den T (e) t dfx N e {E) ( C(E) - dE' ^(E',E) 



dE dt ' J em . J_i 2 \ J E dE 




1 1 .. rE 



den T (e)l d» ^Ji [ dEj ^(^, E) 



1 2 J Epmm ' dE 



Ec - ma * . doc, 

dE 



x[n 7 , in (E) + c dE' N e (E')—^(E',E' - ) (21) 



This equation for electron transport is a restatement of the condition of steady state equilibrium. 

In (f2Tj) AT e is present on both sides. Also the dependence on N e is incorporated in several integral terms. The most 
effective strategy for solving (f2Tj) is to apply the Runge-Kutta technique. 

VII. NUMERICAL TECHNIQUES 

The Runge-Kutta technique is a method for iteratively solving a differential equation which has the following form: 

Vn+l = Vn + hf(x n ,y n ) 

f( x y)-^ 

OjXji 

The solution is advanced through small steps h, with x n +i = x n + h, and the result is accumulated. An adaptive 
step size is incorporated to reduce computation time while holding relative errors fixed. The C++ code utilizes 4 
order Runge-Kutta method with 5 th order error checking. It is prudent to recall that higher-order is not synonymous 
with either smaller error or improved numerical stability. However this technique converges rapidly for (|2"Tj) . The 
details of the Runge-Kutta calculation are described in many texts including Numerical Recipes \3a} . 

I hold absolute errors to 1.0 x 10~ 4 at machine precision and relative errors to 0.0. The equation for electron 
transport gives solutions for the expected number of particles in a logarithmic energy bin. The initial electron 
spectrum is null and the initial 7-ray spectrum is set to unit height in a logarithmic energy bin corresponding to the 
injection energy. The energy range is divided into an arbitrary number of intervals of constant logarithmic width and 
(|2ip is repeatedly solved starting at higher energies and moving to lower. 

The highest energy bin must be solved first since pair production occurs in the highest bins first, and then subsequent 
electrons downscatter to lower energies over repeated iteration of the cascade process. In other words the electrons 
first appear in the higher energy bins and move down in energy. 

The result of this iteration is an electron spectrum corresponding to the specified photon injection. One final 
generation of inverse Compton scattering reveals the outgoing 7-ray production. The output production is integrated 
from the inverse Compton scattering rate to give the final result [2l| . 



n 7 (£ 7 ) = f '"° X dE' N e {E')T KN (E', E' - E 7 ) (22) 

The 7-ray production is then integrated over the time variable t until the output is fully saturated, i.e., until the 
observable energy is equal to the injection energy. The solution of (|21|) is depicted in Fig. O Finally, flux may be 
determined from ([5]). 

Applying numerical solution to (j2"Tj) shows that during the cascade height (or equivalently energy) is preserved 
per logarithmic energy interval on an E 2 dN/dE plot as 7-rays cool. Fig. [5] shows that the cascade process recycles 
most of the injected energy until it finally turns off near E c , then the resulting spectrum is simply the portion due to 
inverse Compton scattering. Fig. [5] uses a logarithmic scale where each range represents a fixed interval of constant 
energy. Energy can be read as the height of the figure. There are a total of 50 constant logarithmic intervals, however 
this number is chosen arbitrarily to minimize processing time while accurately representing the shape of the output 
spectrum. 
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FIG. 5: A plot of fl22} with solutions of J2T| as input. The 50 TeV (top left), 100 TeV (top right), 500 TeV (bottom left) and 
1000 TeV (bottom right) 7-rays are injected (shown in green) above Eth- The spectrum resulting from the cascade has two 
components. The observable portion below Eth in black and the continuously propagating portion above Eth is depicted in 
red. Just below the threshold energy of 10 TeV, the height of the resulting spectrum is about 1/2 the height of the injection. 
Eq. ((5} is superimposed on the 50 TeV injection for comparison. 



VIII. RESULTS FOR CASCADING PARTICLES 

As high energy 7-rays are injected above E t h, they immediately form pairs and begin cycling through the cascade 
process, finally just above E c and after many generations, cascading electrons and photons have nearly the same 
energy as the first generation did. Suddenly, the cascade turns off and most of the energy is preserved in the range 
E c to Eth- 

Therefore, the cascade acts as a photon calorimeter. About 85% of the energy injected above E t h is ultimately 
observable in the energy decade (Eth/W, Eth) near E c . However, the cascade does not preserve the conformation 
of the injected spectrum. While a narrow bin may be injected, the observable spectrum is fixed by the shape of the 
inverse Compton scattering rate and is not accurately reflected as a simple energy translation of the input bin. 

At the factor of two level, one may consider a rule of thumb. 



This is a naive heuristic which hides the relevant physics, however the conclusion that the majority of the injected 
power scatters into the energy range at E c is valid. Essentially this relation demonstrates that the cascade conserves 
energy. 

Importantly the energy due to ultra-high energy 7-ray injection processes is not lost in the sinuous cascade. The 
total output power is equal to the total input power (conservation of energy) and the rapid ascent of the inverse 
Compton scattered photon spectrum ensures that most of the observable energy will appear at E c . These results 
clearly rule out a "bin-shifting" approach to 7-ray energy conservation, however they do provide a number of useful 
estimates of expected spectral outcomes. 

Finally, the cascade takes place extremely rapidly on cosmological scales. Direct computation of redshift effects 
on cascade spectra is not needed since from the standpoint of cascading particles the universe is essentially flat and 
static. One may however easily extend the spectrum to produce real cosmological flux limits using ijTj). 

IX. SUMMARY 

In conclustion, 7-ray energy injected by ultra high-energy isotropic cosmological processes is observable as a spec- 
trum of cooled 7-rays. These processes differ from point sources in that point sources suffer attenuation on the 
cosmological medium. 

The consequences of this work may be summarized as follows: 

1. Diffuse isotropic injection processes above 10 TeV are constrained by experiments which observe 7-rays below 
lOTeV. 

2. The total input power spectrum is reprocessed leading to a 7-ray pile up at E c . This output spectrum gives 
total integrated limits on injection energy. 

3. Observations of present day fluxes suffer important restrictions from both energy and redshift scales. Cosmo- 
logical processes that contribute to the EGRET diffuse flux observations must be significant at z w 1. The 
significance of any observation of present day cascade radiation processes is highly suppressed by a factor of 
1000. 

There are energy loss processes which I have not considered in these calculations. First the synchrotron losses 
may have unit order corrections on these calculations where large magnetic fields are present. Next ionization energy 
losses in dense galactic regions may play a key role in 7-ray attenuation. Finally, bremsstrahlung losses may give a 
measurable correction to galactic fluxes. 

For all of these loss processes, the rate of loss could easily be measured by considering the emission of well understood 
sources. In particular, through the synchrotron mechanism this may give a testable method of determination of the 
intergalactic magnetic fields. 

Throughout this discussion I have taken it for granted that the photon injection was continuous and isotropic on 
the cosmological scale. Individual point sources represent a dramatically different type of propagation problem than 
what is covered in this work. 
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APPENDIX A: THE SPECTRUM OF 7-RAYS SCATTERED AT E c 

Berezinksy derives the spectrum of 7-rays at E c in his textbook [12j | . In the high energy limit electrons conserve 
energy, Eq e (E) — const, above E c . The production of inverse Compton 7-rays is the number in a logarithmic interval 
Qe{E c ) dE over the width of a logarithmic interval dE e / E~ p i.e., dlnE. The spectrum is: 

dE e 

ni = q e {E c )-=^ (*) 

But as we have said, energy is conserved, therefore at E c , q e = qo. The energy carried by an outgoing 7-ray is the 
fraction liberated from the electron or positron, E 1 = fE e , by Eq. J3J: 

E - 4 ^1 
7 ~ 3 ml 

Taking the derivative and isolating dE e we have: 

1 3 ml dE. 



dE„ 



7 



2 4 e E e 



And solving for E e : 



3E 1 



Finally we can combine these relations in Q: 



1 3 dE~ 

n 7 = -q ome] j--^ 



As expected, at E c the spectrum of 7-rays is falling as E 3 / 2 . 



APPENDIX B: THE SPECTRUM OF 7-RAYS SCATTERED AT HIGHER ENERGIES 

Above, we present the conventional discussion of 7-ray spectra resulting from a single cascade step. It is trivial to 
extend this discussion to include an additional step. If the injected particle energy is taken as AEth rather than Eth 
an additional scattering step results. Here the low energy approximation to energy loss fraction does not apply. In 
this "middle" energy range the outgoing electron and positron share total energy so the loss fraction is now: 




Then the following expression trivially follows from J*J: 



dE 1 

n 7 = 2q -^ h "7,2 

In this relation n 7j 2 will be the flux produced by the subsequent outgoing 7-ray produced at E t h- Yet we have already 
deduced this spectrum in App. [XJ Therefore, 




4e E* /2 



The E 1 spectrum will dominate for E > 1. 



